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Abstract 

We obtain the corrections in l/n and in l/lnn (n is the principal quantum number of the 
bound state) of the decay constants of scalar and pseudoscalar currents in two and four 
dimensions in the large N^. We obtain them from the operator product expansion provided 
a model for the large n mass spectrum is given. In the two-dimensional case the spectrum is 
known and the corrections obtained in this paper are model independent. We confirm these 
results by confronting them with the numerical solution of the 't Hooft model. We also con- 
sider a model at finite and obtain the associated decay constants that are consistent with 
perturbation theory. This example shows that that the inclusion of perturbative corrections, 
or finite effects, to the OPE does not constrain the slope of the Regge trajectories, which 
remain a free parameter for each different channel. 
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1 Introduction 



Since the pioneering work of Ref. [Ij, the operator product expansion (OPE) has been 
intensively used in order to improve our understanding of the non-perturbative dynamics of 
the hadronic spectrum and decays for large excitations [2l[3lll[5l[6l[71[8l[9l[l0l[IIl[l2l[T3l[Tl]. 
Specially intense has been the study of the vacuum polarization correlator with different local 
currents. Typically, the large N^. limit [15] is considered in these analysis, since it allows one 
to have infinitely narrow resonances at arbitrarily large energies. 

Recently, this problem has been revisited in Ref. |16j. The aim of that paper was to 
go beyond previous analysis by means of the systematic incorporation of the perturbative 
corrections to the lowest order (parton) result, as well as the subleading terms in 
in the OPE expression of the current-current correlator. In order to reproduce the OPE 
expression from the hadronic version of the current- current correlator, it was necessary to 
include corrections in 1/n (n is the principal quantum number of the bound state) in the 
expressions of the mass spectrum and decay constants. By doing so, it was shown that (in the 
large Nc limit) the combination of the OPE with the knowledge of the large n mass spectrum 
allows one to perform a systematic determination of the preasymptotic corrections in 1/n of 
the decay constants. Even more so, power-like 1/n corrections can only be incorporated in 
the masses and decay constants after the inclusion of the perturbative corrections in as in 
the OPE expression. 

One of the aims of this paper is to continue this line of research. We first consider the 
scalar/pseudo-scalar correlators (in Ref. [T6j only the vector and axial- vector correlators 
were considered). Besides performing the phenomenological analysis, we will have to deal 
with non-trivial anomalous dimensions, which make the analysis significantly different. We 
then consider the 't Hooft model [17] (the large Nc limit of two-dimensional QCD). We 
compute for the first time the 1/n corrections to the decay constants in this model. We 
find that our results agree quite well with those obtained from a numerical evaluation using 
the 't Hooft equation. This provides us with a non-trivial check of our analytical result. 
Finally, we also consider a model at finite N^. and obtain the associated decay constants that 
are consistent with perturbation theory. This model is inspired in the large Nc limit, with 
small widths for low n resonances, but still with the right analytic properties in the complex 
plane. The construction of such a model is not completely trivial. To our knowledge it has 
first been used in Ref. [2]. Here we find that this model survives the test of the inclusion of 
the perturbative logarithms of Q^, though the decay constants are no longer trivially related 
with the perturbative expression of the imaginary part of the correlator. We also show that 
the inclusion of perturbative corrections, or finite Nc effects, to the OPE does not constrain 
the slope of the Regge trajectories, which remain a free parameter for each different channel. 

As we have mentioned above, in this paper we obtain the 1/n corrections to the decay 
constants in the 't Hooft model. This computation is actually one of the main results of this 
paper. It is model independent, since the mass spectrum is known. These corrections are 
an unavoidable ingredient to undertake the analysis of the preasympotic corrections to deep 
inelastic scattering or B decays in two dimensions. Those processes can be in principle be 
studied through the operator product expansion. In order to search for duality violations 
it is compulsory to know the values of the decay constants for large n. We expect to profit 
from our results in the near future to improve over the analysis of Ref. [TB], as well as to 
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perform a similar analysis for deep inelastic scattering [19] . 



2 Scalar/pseudoscalar correlator in D=4 dimensions 

The scalar/pseudoscalar correlator reads 

n5,p(g;/i) = i j d^xe'^^'ivaclT | J5,p(a;)4,p(0)} , (1) 

where Js = "^fi^j and Jp = iipfj^ipf- The dependence on the renormalization scale /x 
appears because the currents have anomalous dimensions. In order to shorten the notation, 
we will use X = S, P in what follows when the distinction is not important. Using dispersion 
relations the correlator can be written like 

f°° 11 

nx(Q';/i)=/ dt——-lmUx{t;iJ,), (2) 

where = —q^ is the Euclidean momentum. In four dimensions the parton-model result 
gives 

nx(g;/i)~g'ln% . (3) 
In order to avoid these divergences, one typically considers derivatives of the correlators like 
. N ^2 d nx(Q^/i) -nx(0;u) , 1 1, ^ , , 

or 

Bx{Q ; /i) = -y (^g2)2 = Q X (ITW- ' ■ 

Since we are working in the large limit, the spectrum consists of infinitely narrow reso- 
nances with mass Mx(n), and these functions can be written in the following way 



where Fx = {vac\Jx{0)\n), and 



n=0 

For definiteness, in this article we will work with the function Bx{Q^', f^)- For large positive 
Q^, one may try to approximate this function by its OPE, which in the chiral and Nc ^ oo 
limit has the following structure 



Bx,oPEiQ'; /i) = ( 1 - \n^^^ (f^ ' + 
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167r2 



:i + (^(«s)) 



:i + (3{a,{i^)){vac\G^{fi)\vac) + 0{^ 



The term /3(G^) above is renormalization group invariant and (3 = ^'^da/ {dix^). 
In order to obtain the perturbative piece of this expression one can use 

imnr*-(t;/i) = ^iv,i?(t;/i) , (9) 

where the expression of R{s) = 1 + 0{as) can be found in Ref. [201 EI] with four loop 
accuracy. The one-loop expression for the coefficient of the gluon condensate can be read 
from Ref. [22]. 

It should be stressed that i?(s; /x) is renormalization-scale dependent and it fulfills the 
following equation (it is understood that i?(s; /i) is also computed in the MS scheme) 

^^'^^ml^{^^)R{s■,^^) = Q . (10) 

Therefore, Imnx(t;/i) does not have a physical meaning by itself (in physical processes it 
should appear associated to masses or equivalent). 



2.1 Matching 

High excitations of the QCD spectrum are believed to satisfy linear Regge trajectories: 

lim — = constant. 

n^oo n 

For generic current- current correlators, such behaviour is consistent with perturbation 
theory in the Euclidean region at leading order in if the decay constants are taken to be 
"constants" , ie. independent of the principal quantum number n. 

In this paper, we would like to include power-like corrections in ctg and 1/Q^ in a system- 
atic way. In order to do so we will have to consider corrections to the linear Regge trajectories 
as well as to the decay constants. We will follow the same procedure used in Ref. [16j for 
the vector and axial-vector channels. We will consider that the large n expression for the 
mass spectrum can be organized within a 1/n expansion in a systematic way starting from 
the asymptotic linear Regge behaviour. In order to fix (and simplify) the problem we will 
assume that no Inn term appears in the mass spectrunj^. Therefore, we write the mass 
spectrum in the following way (for large n) 

MUn) = J2 B^r^n^-'^ = B^xn + bP + + ■■■ (11) 

s=— 1 ^ 

where -B^ are constants. We define Mx Loi''^) — B^n, Mx Nioi''^) — B^n + and so 
on for the leading order (LO), next-to-leading order (NLO), etc. To shorten the notation, 
we will denote B^ = Bx, Bf = Ax and B^^^ = Cx- 

^This is a simplification. If one considers, for instance, the 't Hooft model [TJ, Inn terms do indeed 
appear, as we will see in the next section. If we relax this condition one can only fix the ratio between the 
decay constant and the derivative of the mass. Actually, this can be done in a model independent way. The 
explicit formulas are shown in the appendix. 
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For the decay constants, we will have a double expansion in 1/n and 1/ Inn. 



The logarithmic dependence on n of the coefficients Fx s{n; fi) has the following typical 
structure (70 is defined in sec. 12. 2p : 



/ 1 \ 270 00 

\ / ,, — n 



=0 

As we did with the masses, we will define F|- j^q (n; /i) = F|- g (n; /x) , F|- jy^Q ('^; /i) = F]^ q (n; /i) + 
F|- ;^(n; ii)/n, and so on. Note that in this case we also have an expansion in 1/ Inn. 

We are now in the position to start the computation. Our aim is to compare the hadronic 
and OPE expressions of Bx{Q^] /i) within an expansion in l/Q^, but keeping the logarithms 
of Q. In order to do so we have to arrange the hadronic expression appropiately. Our 
strategy is to split the sum over hadronic resonances into two pieces, for n above or below 
some arbitrary but formally large n* such that Aqcd ^ Aqcd^^* -C Q. The sum up to n* 
can be analytically expanded in and will not generate In terms (neither a constant 
term at leading order in l/Q^). For the sum from n* up to 00, we can use Eqs. f llll) and 
f[T^ and the Euler-MacLaurin formula to transform the sum in an integral plus corrections 
in l/Q"^. Whereas the latter do not produce logarithms, the integral does. These logarithms 
of Q are generated by the large n behaviour of the bound states and the introduction of 
powers of 1/n is equivalent (once introduced in the integral representation, and for large n) 
to the introduction of (logarithmically modulated) corrections in the OPE expression. 

Therefore, by using the Euler-MacLaurin formula, we write i3x(Q^;/i) in the following 
way {B2 = 1/6, B4 = -1/30, ...) 



""^^ Fl(n;/x) r* FUn;fi) 

2-. [Q^ + Mj,{n)Y io (g^ + Mi(n))3 

(14) 



.Q' Fl{n*-fi) - ^(2^-1) Fl(n;^) 

2 (g2 + Mi(n*))3 ' (2A;)!c/n(2fe-i) (Q2 + ]Vf2(n))3 

where n* stands for the subtraction point we mentioned above, such that for n larger than 
n* one can use the asymptotic expressions fllip and fll2p . This allows us to eliminate terms 
that vanish when n — 00. Note that the last sum in Eq. (I14p is an asymptotic series, and 
in this sense the equality should be understood. 

Note also that for n below n*, we will not distinguish between LO, NLO, etc... in masses 
or decay constants, since for those states we will not assume that one can do an expansion 
in 1/n and use Eqs. ffTT]) and ffT^ . 

Finally, note that the expressions we have for the masses and decay constants become 
more and more infrared singular as we go to higher and higher orders in the 1/n expansion. 
This is not a problem, since we always cut off the integral for n smaller than n*. Either 
way, the major problems in the correlator would come from the decay constants, since, in 
the case of the mass, effectively acts as an infrared regulator. 
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And a final comment on renormalons: in principle, we are now in the position to match 
Eq. f[T^ with Eq. ([S]) order by order in l/Q"^. One should keep in mind however that each 
order in of Eq. ([H]) suffers from renormalon ambiguities. Those ambiguities cancel 

between different orders in the expansion. A way to deal with this problem is to devise 
a scheme of subtracting renormalons from the perturbative series, passing the renormalon 
to the condensates, or to higher order terms in the l/Q^ expansion, where the renormalon 
ambiguity cancels (see for instance [23] for an example of such a renormalon subtraction 
scheme). We will not do this explicitly in this paper, since it goes beyond our purposes and, 
with the precision we aim at here, these effects do not appear to be numerically dominant, at 
least for large n. Nevertheless, it remains to be seen whether they lead to some improvements 
for low n. 



2.2 LO Matching 



We want to match the hadronic, Eq. f|T4|) . and OPE, Eq. ([8]), expressions for Bx{Q^',fJ') 
at the lowest order in l/Q"^. Only the first term in Eq. f|T^ can generate logarithms or 
constant terms that are not suppressed by powers of so this is the term that has to 

be matched to the perturbative part of the OPE expansion. We have to consider the lowest 
order expressions in 1/n for the masses and decay constants, i.e. F]. j^Q^n; /i) and M|^^q(?t,), 
since the corrections in 1/n give contributions suppressed by powers of The matching 

condition is then 



dn 



(Q2 + Ml^oin))^ 167r2 



11-6 In 



which can be fulfilled by demanding 



N, 



TT 



+ . . . 



(15) 



\dM'j^ ^Q{n) / dn\ n 



imn^-(M2 H;/^) 



(16) 



This leads us to the following expression for Fp^^ g(n; fi) with four- loop running precision 



B 



X 



87r2 



a{nBx)Y' c{a{nBx) 



c(a(/i2)) _ 



(17) 



(1 + ria{nBx) + r2a{nBxf + na{nBxf) 
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where a(/x^) = as{fj,'^)/n, 



a{nBx) = a(/i^) s 1 + a{^^)j3Q In 



nB 



X 



a(/i2)3(/5oln(^ 



+ In (1 + a(/i2)/5oln 



1 + a(/i2)/5o In ( ^ j + a(/i2)/5i In M + a(/i2)/3o In 



+a(/i^)/3iln 



1 + a(/i^)/3o In 



/ l + a(/i2)/5oln(^f^ 
+a(/i2)/fi In ( 1 + a{fi^)(3o In ( ^ ) + a{fi^)Pi In ( 1 + a{^^)f3o In 



2a(//2)4/?oln ( nfc) + a(/i2)5/52in2 (iigK 



l + a(/i2)/5oln(^ 



and 



c(x) = 1 + (71 - /3i7o)a; + ^ (^(71 - /5i7o)^ + 72 + - Pili - /527o) a;^ 
+ Q(7i - Pilof + ^(71 - A7o)(72 + A^To - Pill - 



+|(73 - ASo + 2/3i/327o - f^slo + I3ili - I32li - A72) ) x"^ 



^ 2 



(19) 



with 7j = 7j//9o and /3j = Pi/Po- The values of the various constants are, in the large N^. 
limit, 



ri 

^^3 



l7Nr. 



/V2 

;^(7431 - IGOtt^ - 1920C(3)) 



(20) 



497664 



(25999999 - 10952647r2 - 11200032C(3) + 1607040C(5)) 



iijv n 17A^c a _ 2857Ar3 ^ (150553 - 2376C(3)) 



^° 12^ ' /^i 24' ' 3456 ' ' 124416 

185N^ 11413Ar3 A/'^4(460151 + 74048C(3) - 126720C(5)) 

7o = -g^ , 7i = , 72 = ,000. ' 73 



384 ' ' 13824 ' ' 294912 



The MS mass anomalous dimensions are taken from Ref. [211 [25] . 

Finally, we remind that, strictly speaking, we can only fix the ratio between the decay 
constant and the derivative of the mass. We have fixed this ambiguity by arbitrarily imposing 
the n dependence of the mass spectrum. 
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2.3 NLO matching 

We now want to obtain extra information on the decay constant by demanding the vahdity 
of the OPE at in particular the absence of condensates at this order. We insert 

the NLO expressions for Mx{n) and Fx{n] fi) into Eq. f|T^ and impose that there be no 
l/Q"^ contribution. With the ansatz we are using for the mass at NLO, it is compulsory to 
introduce the (logarithmically modulated) 1/n corrections to the decay constant if we want 
this constraint to hold. Note that it is possible to shift all the perturbative corrections to 
the decay constant. 

Imposing that the term vanishes produces the following sum rule: 



1 ^, ^.MA 



n*-l „„* 
n=0 "^0 



oo 



dn 



2Q2 



(g2 + M2 (n))3 



+Q / dn ^^^ , 'V, = . (21) 



(Q2 + Ml,o{n)) 

This equality should hold independently of the value of n*, which formally should be taken 
large enough so that as{n*By) <^ 1, i.e. the limit <^ n*By <C Q^. Only a few terms in 
Eq. (I2TI) can generate In terms, which should cancel at any order. Those are the first two 
and the last two terms. Actually, the next to last term does not generate logarithms, but 
it allows to regulate possible infrared divergences appearing in the calculation. Therefore, 
asking for the cancellation of the suppressed logarithmic terms produced by the first 

two and the last term in Eq. (12T!) fixes Fx 1 . The non- logarithmic terms should also be 
cancelled but they cannot be fixed from perturbation theory. 

One can actually find an explicit solution to the above constraint for Fy^ by performing 
some integration by parts. We obtain 

Fh{n;^^) = ^^{nFl,{n;f,)) . (22) 

Note that Fxi{n]fi) is of the same order in as as F|-Q(n;/i). This is different that in the 
vector/axial- vector case. For further details to the procedure we refer to Ref. [16]. Here we 
would only like to insist that non-logarithmic terms should also be cancelled, but they cannot 
be fixed from perturbation theory. For these terms we cannot even give a closed expression: 
In Q^-independent terms may receive contributions from any subleading order in the 1/n 
expansion of the masses and decay constants. The reason is that the decay constant at a 
given order in 1/n is obtained after performing some integration by parts, which generates 
new (In Q^-independent) terms that can be enhanced. This statement is general and also 
applies to any subleading power in the 1/Q^ matching computation. 

2.4 NNLO matching 

We now consider expressions for the mass and decay constants at NNLO. For the first time 
we have to consider condensates. Simplifying terms that do not produce logs, we obtain the 
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following equation, 

1 3 /as(Q'^^"'° 



iV,Q2 22 Vas(/i2 



P{asiiJ,)){vac\G\lJ,)\vac) (23) 



where = stands for the fact that this equality is only true at leading logarithmic order. Using 

270 POO -1 1 / ^ f ^^\270 



we have 



n4 1 K- f as{Bxn) yPo 



Fhin;^) = ^^(^lA^Fl,in;fx) + CxFl,{n;f^)^ (25) 
3 _ /'as(5xn)y'^° /5(as(/i))(^^ac|G'2(/i)|t;ac) 

Note that the accuracy of this result is set by our knowledge of the matching coefficient of 
the gluon condensate. Note as well that Fjr^{n] jj) is suppressed with respect to F|- ^{ri] jj) 
and F^o{n;fi). 

2.5 Scalar versus pseudoscalar correlators 

From the previous analysis it is evident that the coefficients of the mass spectrum are free 
parameters and can be different for the scalar and pseudoscalar channel, in other words, 
they cannot be fixed from the OPE alone. This point was already emphasized in Ref. jH] 
for a model that reproduces the parton-model logarithm. In this paper, we show that that 
the inclusion of corrections in as does not affect that conclusion, and that Bs and Bp 
somewhat play the role of the renormalization scale in the analogous perturbative analysis 
in the Euclidean regime, and are therefore unobservable. Overall, the situation is similar to 
the case of vector and axial- vector correlators studied in [3 [16]. 

However, although the constants that characterize the spectrum can be different for the 
scalar and pseudoscalar channels, they have to produce the same expressions for the OPE 
when combined with the decay constants. This produces some relations that must be fulfilled. 
Defining t = Bsn = Bpn' and taking n and n' as continuous variables, these relations are 

^^Fl,{n■^^) = J^fo{t■^^) , (27) 

n d f\ \ 
Fl2in;fi) - — ^ i^-AxFl,in;fi) + CxFl,in;fi)j (28) 

3 _ / a,(t) ..p{as{i^)){vac\G\iJ,)\vac) 



-^Fl,{n; /i) = --Imnf ^•*-(t; /i) ^ ^fo{t; /i) , (26) 
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2.6 Numerical analysis 

The aim of this section is not to perform an in-depth numerical analysis of the expressions 
we have found. We lack experimental data for resonances at high excitations, where our 
expansion in 1/n would work best, and have no information on the decay constants; also, 
one must not forget that we are staying at leading order in the 1/Nc expansion and in the 
chiral limit. Instead, the purpose of our analysis is just to get a feeling of the relative size of 
the corrections in l/ra, and of the importance of the resummation of the powers of in f[T7|) . 
We restrict ourselves to the SU(2) flavour case. In the table we show, in parenthesis, the 
experimental values of the masses of the mesons. We fit these values to Eq. (ITT]) , and from 
these fits we find the decay constants. We take most of the experimental values from [27J, 
except the pseudoscalar states with n = 3,4, which we take from [28]. We do not include 
the pion (the state with = 0) as input in our fit. We define 

F|p(n; /i) = 2Bs,pnGlp{n- /i) . (29) 





n = 1 


n = 2 


n = 3 


n = 4 




986(980 ± 10) 


1342(1370) 


1544(1507 ±5) 


1703(1718 ±6) 




1305(1300 ± 100) 


1791(1812 ± 14) 


2098(2070 ± 35) 


2349(2360 ± 30) 


Gs 


13318 


452 


229 


188 


Gp 


649 


301 


255 


233 



Table 1: We give the experimental values of the masses (in MeV) for scalar and pseu- 
doscalar particles (within parenthesis), compared with the values obtained from the fit. We 
take as{l GeV) = 0.5 and p{G^) = -(352 MeV)^ Note that the values ofGs and Gp depend 
on the factorization scale. We have taken fi^ = lOBs and fi'^ = lOBp for the scalars and 
pseudoscalars respectively. 

The parameters of the mass spectrum obtained from the fit to the experimental values 
in the table are 

^5 = 0.456 GeV^ As = 1.262 GeV^ = -0.746 GeV^ (30) 

5p = 1.040GeV^ Ap = 1.589 GeV^ Cp = -0.926 GeV^ (31) 

In Figured] we can see the difference between the resummed expression of G'j^ i^Q{n; fi) , 
as taken from Eq. ( flTl) . and the expanded expressions, at different orders in ag. We can see 
important differences with respect to the analysis of Ref. for the vector and axial- vector 
currents. The impact of the resummation of logarithms (or of perturbation theory in general) 
is much more important for the scalar and pseudoscalar currents, in particular for the former. 
Note as well that the importance of these perturbative corrections enforces us to take a large 
value for the factorization scale to ensure nice convergence properties of the perturbative 
series. By inspection of the plots we conclude that our figures are a reliable approximation 
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for, at least, n = 2 (being on the boundary) or larger. This should be compared with the 
vector and axial-vector case, in which even for n = 1 one obtains reasonable numbers. One 
should keep in mind, however, that the decay constants we are computing here are scheme 
and scale dependent quantities. Therefore, it would be more appropriated to consider them in 
combination with another quantity with the inverse scheme and scale dependence to become 
a direct observable quantity. 

Figure [2] shows us that the dependence of the decay constants on n is small (again from 
n ~ 2 on). The corrections in 1/n are rapidly converging, the result at NNLO being almost 
indistinguishable from that of NLO in the region where the series is convergent. One may 
ask whether these results change qualitatively by choosing a different value of the gluon 
condensate. We remind that the value for the gluon condensate will vary depending on the 
renormalon subtraction scheme used. We have seen that the change is quite small if we vary 
the gluon condensate between the range 0.01 and 0.04 GeV"^. 



Gs,Lo(n;10Bs) 
350- 
300- 
250- 
200- 
150- 
100- 
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Gp,Lo (n; lOBp) 
325- 
300- 
275- 
250- 
225- 
200- 
175_ 
150- 




Figure 1: In this plot we show the differences between the expanded and the resummed ex- 
pression for G(^s,P),Lo{^i yu); following the definition given in eq. ^EW- 




11 



3 Preasymptotic effects in 1/n in the 't Hooft Model 



A systematic analysis of the preasymptotic effects in 1/n has not been undertaken until the 
analysis of Ref. |16j and this paper. Therefore, it is interesting to check the methodology in a 
specific model where the results obtained with the OPE can be tested with the exact results. 
Here we consider the 't Hooft model, which will allow us to perform such an analysis. In this 
model, the preasymptotic corrections to the mass spectrum are known [17] (actually they 
have a In n dependence, which will prove to be crucial) but not those to the decay constants, 
for which the analysis of the preasymptotic effects in 1/n is not an easy issue. 

Unlike in the rest of the paper, we will consider the general situation with non-zero (and 
different) quark masses. We will focus our interest in the scalar and pseudoscalar currents. 
One reason is that in two dimensions the vector and axial- vector current matrix elements can 
be related to them. Nevertheless, the main motivation is that we are interested in situations 
where logarithms of Q can be generated, which can only happen if an infinite number of 
resonances contributes to the correlator. This does not happen for the vector and axial- 
vector correlator in the massless limit, where they both become almost trivial, since only 
one resonance (the ground state) has overlap with the currentl^ P^ IBU] . 

3.1 Preasymptotic effects in 1/n from the OPE 

The main aim here is to repeat up to NLO and in two dimensions the analysis performed in 
sec. 121 The difference is that now we know the spectrum at NLO, which reads [T7] 

oo 

Ml{n) = M\n) = ^ B^-'^n^-'^ = B^^^n + 5^°) + ■ ■ • , (32) 

and we define MlQ{n) = B^^^n. Note that the spectrum is the same for the scalar and 
pseudoscalar channel. The coefficients B^^^ and B^^^ are known in the 't Hooft model. They 
read [ITj 

S(^) = TT^P^ , 5^°) = (mj^ + m]^^) Inn + constant , (33) 

where i and j represent, respectively, the flavor of the quark and antiquark that make up 
the meson in the 't Hooft model, = g'^Nc/{2TT), and nif^ = nif — is the renormalized 
mass. The explicit expression for the constant term can be found in Ref. [32] . 
Similarly to the mass we write 

Ft 1 (n) 

Fl(n) = Fl.o(n) + ^^ + ..., (34) 

n 

and Fjr ^Q{n) = F^Q^n), Fjrj^^Qiji) = F^Q{n) + F'^ ^{n)/n. Note that, unlike our former 
expression for the decay constants, Eq. f[T^ . there is no n multiplying in front. This is due 
to the fact that in two dimensions the correlator Ux goes like 

Hx(g)~ln^. (35) 

^In fact, the combination of these results with the behaviour of the ground state wave function in the 
X — > hmit (or in other words its mass) is used to fix the quark condensate |31) . 
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The correlator has also no anomalous dimension in the 't Hooft model. The definition of 
the correlators used in this section is the same to the one given in Eq. ([T]) but with general 
currents (with in principle different fiavours): Jg = ipjipi and Jp = ii^j'-j^ipi. The hadronic 
expressions for the correlators read 

oo _2/ \ 

n=l,3,5... ^ ' ^ 



r,—r\ OA \ ' ^ 



where Fx{n) = {vac\Jx\ij]n). We note that only odd and even states contribute to the 
scalar and pseudoscalar correlator respectively (this implies that for the scalar correlator 
the ground state does not contribute). This result was obtained in Ref. [29j, using some 
general properties of the 't Hooft ecmation operator. Those properties can be understood 
somewhat on a pure symmetry basiqj. Nevertheless, it is not an easy task to visualize this 
symmetry at the Lagrangian level. The situation is similar to the one found in the spectrum 
of diatomic atoms [33] or of NRQCD in the static limit (see for instance |S1)- In those 
situations a convenient way to deal with the problem is to project the Lagrangian to the 
Hilbert space sector we are interested in. In our case, this is basically equivalent to ending 
up in the 't Hooft equation (wich we present in the next section), where one can use the 
operator identities found in Ref. [2^ for the current matrix elements (however it should 
also be possible to relate these operator identities with some underlying symmetries of the 
system), in particular Eq. (|50|) . Using these results one can discriminate which states give 
non-zero contribution to each correlator. 

We now define the Adler-like function in two dimensions 

. . (38) 



for which its hadronic expressions read 



nx 



where nx stands for the sum for the scalar or pseudoscalar case respectively. 

In order to follow the procedure used in the previous section we need the OPE expression 
for Ax- It reads: 

AT- = -(l + ^ + Dx^(\n(^)+ln(^))+ ^z,,-^(^^^^) + ^^^^^^^^ 



+ ■■■, (40) 



■^We cannot use parity symmetry due to the fact that the quantization frame that we use does not 
respect this symmetry exphcitly. This means that the states and currents have complicated transformation 
properties under this discrete symmetry. 
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Figure 3: Diagrams that contribute to the OPE at one loop in perturbation theory. 

where Ds = 1 and Dp = —1, and we have neglected terms of 0{m'^ / Q"^) (except for the 
logarithm term) and terms of 0{1/Q'^). The pure perturbative piece of this result is obtained 
from the computation of the diagrams shown in Fig. [31 The condensate contribution is 
obtained by "collapsing" one of the quark propagators. We see that the coefficient of the 
quark condensate has no InQ piece. 

Note that in two dimensions the coupling constant has dimensions. Therefore, the 
perturbative and OPE expansions are related. In particular, in two dimensions, there are 
no problems with renormalons, due to the fact that there are no marginal operators in the 
Lagrangian. 

In Eq. fllUl) the quark condensate, as usual, only includes the pure non-perturbative 
contribution, ie. the free result^ 



_/\r 777,2 

{'ipii'i ) pert. = 7"^ mi In ^ (41) 



has been subtracted and displayed explicitly in J^-^^^' . 

We can now combine all this information with our knowledge of the large n behaviour of 
the mass spectrum to obtain the 1/n preasymptotic effects of the decay constants. 

We follow the same procedure of sec. 12. 1[ To rewrite A^'^'^' using the Euler-Maclaurin 
formula, we first define n = 2n' + 1 for the scalar case and n = 2n' for the pseudoscalar case, 
so that we can transform the sums into an integral over a continuous variable. The only 
term that will give a LO contribution, not suppressed by powers of Q^, will be 

io 2 (M2^(n) + Q2)2 2 (n7r2/32 + Q2)2 2n ' ^ ^ 

where the last equality is the matching to the partonic result and gives us 

Fl,{n) = N,n(3' . (43) 

In order to go beyond the parton result we need to consider the 1/n corrections to the 
linear behaviour in the mass. Note that the form of the spectrum is different from the 
one assumed in Ref. [TH], or in sec. [D in the 't Hooft model we know that B'^^^ is not a 
constant and diverges logarithmically in n. We actually do not care about the finite piece 



''Note that we have a sign difference with respect to the result of Ref. |35j . We also disregard finite pieces, 
which are scheme dependent. 
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of B^^\ since such contribution would constrain the term of the decay constant 

but not the 0{l/n) term; this is so because F|-q is constant. Therefore in what follows we 
will neglect the constant term in B^^^ and only consider the logarithmic correction. Within 
this approximation, the matching at NLO reads (at logarithmic order, after expanding and 
integrating by parts). 



dn 



FlAn)/n 



2 (n7r2/?2 + g2)2 



dn 



-^1,0 



NrD 



X 



27r g2 



(2m,m,) ln(Q2) = _q2 



2 (n7r2/?2 + g2)2 

°° dn 1 

* T(n7r2/52 + g2)2 



nir'^f]'^ 



(44) 



^Dx{2mimj) 



where again = means that the equality is only true at logarithmic order. And so. 



n 



„2m,m,i5v 



m 



i,R 



m 



■j,R 



(45) 



It will prove useful in the next section to combine the scalar and pseudoscalar result for 
the decay constant in a single function. We define 



[n 



Fs{n) forn = odd 
Fp{n) forn = even 



Therefore, we have 



[n] 



1 + 



1^ ^ 2mimj ^ 



-I) 



(46) 



Let us note that there is a term which is sign-alternating. 

Finally, we stress that we have obtained these expressions for the decay constant using 
the OPE, symmetries and the knowledge of the mass spectrum, nothing else. In the next 
subsection we check whether the explicit numeric computation of the decay constants from 
the 't Hooft equation fulfills these expectations. 



3.2 P r easy mpt otic effects in 1/n from the hadronic solution 

In the 't Hooft model it is possible to write the decay constants F^{n) in terms of the 
light-cone distribution amplitude of the bound state, 4>l{{x), which is the solution of the 
equation 

+ Yz^j C(^) - I ^y^n{y)^^^^2 ' (47) 

where x = p'^ /P^, with p+ being the momentum of the quark i, and P stands for Cauchy's 
Principal Value. 

The decay constant for the scalar case then reads [29] 

F5 77 = \ dx<p'Hx) — ^— =m,\ — / dx-^^ for 77 odd, 48 

V 47r Jo \ X 1-xJ V vr Jo ^ 
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and zero otherwise. For the pseudoscalar we have 
"1 



Fp(n) 




rrij rrii 
X 1 



X 




TC 







1 (7)«i(x) 
(ix " for n even , (49) 

X 



and zero otherwise. Note that in the above two equahties we have used the remarkable 
identity 



rrii 







X 



dx- 



X 



(50) 



obtained in Ref. [29] . 

By comparing Eqs. (HHl) and (H9|) with Eq. (H6|) we can obtain expressions for dx '^"J^-' 
with 1/n accuracjl^: 







vr- 



/5 



rrii 



1 + 



-ir + o{ \ 



By using symmetries and the 't Hooft equation, we can also obtain 



M^{n) / dx^'^ix) 



dx ' + 



X 







1 + 



m: 



+ m 



1 — X 
miiTLj 



(51) 



(52) 



2n7r2/32 



+ 



+ (-l)"7r/5mj- 



mirrij 



2n7r2/32 n7r2/?2^ 

where in the last equality we have made use of the symmetry property of the 't Hooft wave 
functions, 

0^^'(^) = (-l)>^Xl-^) • (53) 
These results can also be used to deepen our analytic understanding of the 't Hooft wave 
function in the end-point regions in the massless limit (actually we need only one mass to 
go to zero at each boundary, rrii for x — >■ 0, or mj for x — 1). In the massless quark limit 
the integral in Eq. (15T!) is dominated by the behaviour of the wave-function in the x — >^ 
boundary: 

0jf(x) =c„x'^»(l + o(x)) , (54) 

where Pi is the solution of 

mj^ + /5^7r/?j cot 7r/?i = , (55) 
which in the massless limit approximates to 



Pi = 3" + o{mi} 

71 P 



(56) 



^ In this result the global phase has been fixed to 1. One always has this freedom but note that this also 
fixes the global phase of (x) , which is no longer arbitrary. The leading order result was first obtained in 
Refs. [29l[30] and later confirmed using the boundary layer equation in Ref. [32]. The 1/n correction is new. 
Note that these 1 /n corrections have to be included in any analysis of preasymptotic effects in the 't Hooft 
model, in particular in the analysis of Ref. [18]. Nevertheless, the results for the moments obtained in that 
paper remain unchanged. 
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Therefore we obtain for the integral 



TT/3 



hm TJii / dx— = Cn—p^ , (57) 

m,^0 Jq X a/3 

and, in the massless quark hmit, the problem of getting the decay constants could be refor- 
mulated as that of obtaining the coefficient c„. One can find its value for n = (note that 
this result does also require mj = 0): 

lim Co = 1. (58) 

In principle, there are several ways to obtain this result. One can work along the lines of 
Ref. |36] to obtain an approximated Schrodinger-like equation, which can be approximately 
solved for the ground state. Another possibility to fix the value Cq is by matching the solution 
00 = 1 (the exact solution in the strict massless limit) and the solution (pQ = cqX^' in the 
region of overlap (the latter is valid for x <^ 1, whereas it can be approximated to a constant, 
Co, for values larger than e"^/*"*, which is a very small quantity for small masses. Therefore, 
there is a region on which the constant solutions: "co" , and "1" , overlap and should be equal 
by continuity). One can also use the value of limm-^o Jq 4^o{x) = 1, to fix cq. 

Nevertheless, in this paper, we are mainly interested in the study of high excitations. 
The combination of Eq. fl57|) and Eq. fH6|) allows us to give the following prediction for Cn'- 

(59) 



(1 1 \ 
n^n 2mi^(3^ J 



How much of all this can be understood by a direct analytic computation? In principle the 
't Hooft equation can only be solved numerically. Nevertheless, for large n, the 't Hooft 
equation can be approximated by the boundary- layer equation [30] : 

0^(0 = ^0'(O - / «'(e')P77^^ ' (60) 



where 



^ Wo {i'-i? 



lim0^^-(e/M2(n)), (61) 



as far as we are not in the x — > 1 limit (one could also write a symmetric equation convenient 
for x ~ 1 region). It is possible to analytically solve the Mellin transform of this equation 
|32j . In particular, one obtains 

di^ = --. I diHi) = -j. (62) 

From these results one easily checks the leading order terms of Eqs. (|5T|) and Eq. fl59l) . Note 
that for the last equation (pnix) ~ c„x^' is valid for x -C 1. In particular, for high excitations, 
the region of validity of this expression is very small, restricted to the region ^ = M^x <C (3"^. 

In order to check the 1/n corrections in Eqs. fIFIl) and dSHD, one should be able to go one 
order beyond the analysis of Ref. [32], which appears to be a formidable task. This would 
go far beyond the aim of this paper. Instead, in the remainder of this section, we will try to 
see whether Eqs. (1511) . (158|) and (159|) can be confirmed by a direct numerical computation. 

In order to perform the numerical analysis we will use two methods: 
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1. One is based on the Brower-Spence-Weis improvement of the Multhopp technique [32]. 
In this method the 't Hooft wave function is decomposed in a basis of sin functions. 
Therefore, the 't Hooft equation becomes a infinity-dimensional matrix, which at the 
numerical level is truncated for a finite number of sines. This method appears to be 
very well suited for very high excitations, and allows us to work with different quark 
masses without problems. Nevertheless, it has some problems for very low masses and 
it does not have the right functional behaviour in the x — > 0, 1 limit. 

2. The other method that we use is based on the decomposition of the 't Hooft wave 
function in a basis of — x)'^^P{x) functions, where P{x) are Jacobi polynomials. 
This method has already been used in Refs. [23, ES]- Unfortunately, we are only able 
to give reliable numbers for even states with equal masses. For those states we can 
perform numerical checks with the numbers given in [38] (we use them as numerical 
check of our implementation of the method, since our interest is on a different regime 
than that of [38]: high excitations and small masses). In any case, this method does 
not appear to work very well for very high excitations, as well as for very small masses 
(though it can reach lower limits than method 1). Moreover, at a certain point the 
introduction of more Jacobi polynomials in the calculation spoils the convergence of 
the result, except for the wave function of the ground state. On the other hand, by 
construction, this method should have the right functional behaviour in the a; — 0, 1 
limit. 

We now first try to check Eq. fl5T]) . In Fig. IHwe show our results with method 2. The 
results from the numerical calculation are presented with some rough estimation of their 
uncertainty, performed by considering the difference between evaluating directly Eq. flFTl) or 
a properly weighted (see tables [2] and [3] for the definition) combination of Eq. fl52]) . We show 
in tables [2] and [3] the differences for two mass values. With this method we cannot go to 
very large values of n, nor consider different quark masses or odd values of n. Nevertheless, 
for the range of values one can consider with this method the agreement is perfect. The 1/n 
and mass dependence can be unambiguously checked very nicely. 



n 


7 frfx'^-^") 

p J X 


M'^(n) 


/ dx(t)n{x) 


2 


3.24 


3.24 


4 


3.20 


3.20 


6 


3.19 


3.18 



Table 2: Results from the numerical calculation using method 2 for rrii = rrij = (3. 

In Fig. [5] we show our results with method 1. In this case we can take much larger 
values of n, although not so small quark masses as in method 2. On the other hand, in 
this case, the use of Eq. ( 1HT]) or the properly weighted combination of Eq. (15^ is basically 
indistinguishable. We can see clearly in this figure that scalar and pseudoscalar states follow 
different curves, as their corrections involve the difference of the masses of the quarks or their 
addition, respectively. The behaviour of the corrections changes from increasing to decreasing 
with n depending on how the sum and the difference of the quark masses compare to 2(3. 
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n 


p J a; 


M'^{n) 


/ dx(t)n{x) 


2 


2.97 


2.96 


4 


3.06 


3.05 


6 


3.10 


3.08 



Table 3: Results from the numerical calculation using method 2 for rrii = mj = 0.1(3. 

We can also see how the numerical results approach the expected analytic expressions as n 
increases. 

We conclude that we have been able to unambiguously and nicely check Eqs. fl5T]) and 
fl52|) numerically. We note that we have been able to visualize the mass dependence and the 
sign alternating terms, and that, in particular, we are reaching a numerical precision below 
1%. 

We now try to check Eq. (jSSD- In Fig- El we show the plot of (j)Q{x)/x^^ for small values 
of X calculated with method 2, and we can see that, indeed, as we approach the massless 
limit Co tends to 1. Method 1 is not efficient for the evaluation of the wave function near 
the origin, for, as we said before, it does not have the right functional behaviour in the 
X — s> 0, 1 limit: it yields a sinusoidal function that oscillates around the right solution, until 
we approach the boundaries, where the curve that the oscillations follow does not die off as 
it should. 

Finally, we try to check Eq. (1591) . In this case we cannot use method 2: although it 
incorporates by definition the right c„x^' behaviour, its numerical precision at high values 
of n is very poor. We are left then with method 1. Numerical precision in this case is not 
so much of an issue, but due to the aforementioned limitations of this method we can only 
hope to get some estimate for c„ by finding the path around which our sinusoidal solutions 
oscillate when x is not too close to 0, and then extrapolate this path to the limit x ^ 0. In 
figure [7] we show our results for n=4 with = 0.001/5, rrij = 2/5: we evaluate 04(x)/x^' at 
points separated a distance Ax = 10~^, starting at x = 10~^ and ending at x = 10~^, and 
we fit them to a quadratic function. Depending on how many of the points that are closer 
to the origin we include in the fit, the resulting curve will be one or another, as shown in the 
plot, somewhat oscillating around the two extreme curves displayed. The possible values of 
the fit at X = are the possible values of C4 that we can find with method 1. In figure [S] 
we show the results up to n = 25 for the cases rrii = 0.001/5, rrij = 2/5 and rrii = 0.001/5, 
rrij = 0.001/5. Overall, this numerical analysis serves as a consistency check of Eq. fl59l) . 

4 Finite Nc 

This far the discussion has been restricted to the large Nc limit. We would like to finish this 
paper by studying models valid at finite A^^^. One model with the right analytic properties for 
the current-current correlator and with the correct limit in the large Nc limit has been built in 
Ref. [2]. This model is also able to reproduce the leading parton- model logarithm at large Q^. 
Working analogously to previous sections, we will consider whether subleading perturbative 
corrections can be incorporated in this model, and see if by including 1/Nc effects there 
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3.3 



3.2- 




- — = ml r) 

2 ' ' ' ' 3 ' ' ' ' 4 ' ' ' ' 5 ' ' ' ' (T " 



Figure 4: In this plot we show our results for Eq. ( f37]) with method 2. The hands reflect the 
difference between using directly Eq. 137]) or a properly weighted combination of Eq. 15^) . 
which we show in tables\^ and\^ 




5 10 15 20 



Figure 5: In this plot we show our results with method 1 for Eq. ( f37]) . 
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1.15 



1.10 



m,=0.1/3 



1.05 



m,=0.05yS 



1.00 



m,=0.005/3 



0.95 



' 0.002 0.004 0.006 0.008 0.010 

Figure 6: In this plot we show our results with method 2 for Eq. ( f5^) with rrii = mj. Due to 
the small values of the quark masses, in the range we are showing the functions are essentially 
constant. 



appears any connection between the slopes of resonances of different parity. For definiteness, 
we will consider the vector- vector correlator but most of the discussion trivially applies to any 
other current-current correlator; in particular the scalar or pseudoscalar correlator, for which 
the only complication would be to take into account the anomalous dimensions. Therefore, 
we take 

(g) = {q^'q" - ^'^V)ny(g) = i j d^xe'^%vac\J^{x)J'(r{'^)\vac) , (63) 

where Jy = order to avoid divergences, we will work with the Adler 

function 

. -Q^^^v&) = [ *^^ilmn.«, . (64) 
We consider the following expression for the vacuum polarization, 

"=o(f) ^"^A2 + M2(n) 

For Fy{n) =constant and My{n) = Byn, we recover the model of Ref. [2]. Here we will 
allow Fy{n) to be n-dependent. The Adler function for this model reads 

A (n^^ (^ (Q^V'^^ .... 



21 



2.0 




0.0006 / \ 0.0008 \ 



.0010 



Figure 7: The dots in this plot correspond to values of (f)i{x) / x^^ sampled with Ax = 10~^ for 
rrii = 0.001/5, rrij = 2(3. We show some of the possible fits to a quadratic function, depending 
on the set of points used, from all those that fall within the interval x G [10~^, 10~^] to only 
those that are within x G [10""^, 10~^]. The possible values of the fits at the origin give us an 
estimated range for C4. 



For simplicity we will set = By in the following. It will also be enough for our purposes 
to keep the analysis of A{Q'^) to the lowest order in the OPE. Therefore, we only consider 

A.,^,,m-A{^^^:^)--Q-'^A[it^^ . (67) 

where and aM{t^) admit an analytic expansion in terms of as(Q^) (computed in the 

MS scheme): 

aAQ') = a.(Q')fl + £i^5^ + *^5(f!)+...) (68) 

\ Ci vr Ci vr^ / 

aM[t) = as{t) 1 H \ — H . (69) 

The coefficients q and di have been computed in Ref. [26] (they can be obtained from each 
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other through dispersion relations): 

ci = ^ (70) 



C2 



^ (^(243 - 176C(3)) - ^(11 - 8C(3))) 



Ar3 

C3 = [346201 - 29047r2 - 324528C(3) + 63360C(5) 

+ ^(-62863 + 5287r2 + 51216C(3) - 5760C(5)) 



{EfQf 



/176 



c 



E^g2 1024 V 3 



- 128C(3) 



and 



3 167r2 
/ 



As before, through the Euler-Maclaurin formula we can see that the matching between the 
perturbative and the hadronic calculations at logarithmic order requires that 



n J \ ttNcJ Jo Bv {z + n 

where we have added the subscript to Fy{n), in agreement with our notation of sec. 12.11 
In Ref. [16], we obtained that 

^-^-aJi + cT^^^) , (72) 



Bv - " V ^ 
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where the upperscript oo stands for the Nc ^ oo hmit. One could consider that the same 
structure should hold for finite in our model and that the decay constant should then be 
proportional to 1 + Ci °^^'-^^"\ However, this is not so due to the specific functionality on 
that we have introduced. For instance, it is interesting to consider the following quantity 



dn- 



n + z)2 



aM{Bvn) = aA^Byz) = a^iQ"^ 



^ta£/3XWW^))+0(i,).(73) 



The first term is what we should obtain to match the perturbative calculation. If we want 
FyQ{n) to be expressed in terms of as{Byn) we have to fine tune it in order to eliminate the 
second extra term. 
If we write 



P2 



By 1- J: 



ci- 



{Bvn) \ 

TT / 



C 



J ttNc B 



(74) 



V 



the following equality has to be satisfied (here, with the precision of the calculation, we can 
replace z ~ Q'^/By) 



dn- 



1 







ln-|^/5I(«^(Q')). 
By 



(75) 



(n + z)^ By 

This implies that ~ a'^{By) Inn. Note that Pq and Pi are the same for a^, and as- 
Therefore, at low orders in a the distinction between different renormalization schemes is 
superfluous. 

We can try to refine our ansatz. If we rewrite the expression for Fy in the following way, 

c 



A, 



B 



V 



1 



_0 

c 



1 + Ci 



c 

ttNc 




aM{Byn^^^N, 



TT 




nNr B 



V 



(76) 



aM{Byn) 



IT 



^-^lnn/?-(a-(i?v^n)) + 0(^; 



1+: 



some logs are reabsorbed in aM{Byn 
easily obtain the coefficient multiplying a^{B 



and SFyQ ~ 



0{aUBy),a!+^{By] 
We find that, at C(af), 



B 



V 



^ 9 

-—PoCi^siBi 
on 



C5F^ 
nNc By 

In* ra). We can 
(77) 



The calculations for the axial- vector case are equal replacing By by Ba (the slope of the 
Regge trajectory in the axial- vector case). Therefore, as in Ref. [T6], we can satisfy the 
constraints of the OPE with independent values for By and Ba- Thus, we conclude that 
the OPE cannot constrain By and Ba to be equal for general models with the right analytic 
properties at finite Nc- Note as well that, in general, at finite Nc 

FUn) 



1. 



\dM^ 



V,LO 



'n)ldn\ 



^-Imn^^*-(M^,^oH)- 



vr 



(7J 



^Nevertheless, in this case it is more difficult to get a closed expression with as+*(i?v') In* n accuracy. 
We do not undertake the effort to try to get it. At this level of precision our expression is good enough and 
also good enough to make our point. 
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4.1 p at finite A^^ 

We now briefly consider what this model tells us about the computation of p(t) = ^lmllv{t)- 
This is a physical observable related with a{e~^e~ — > hadrons). 

We first note that one can approximate the sum over n by an integral (i.e. use the Euler 
McLaurin expansion) if oo and | arg^l < vr for the model of Ref. p]. At this respect 
the modification to this model introduced by us is simply the introduction of logs of n in the 
decay constants, which should not change the analytic properties of the function. This allows 
us to obtain the OPE expression for the vacuum polarization from the hadronic expression 
if I argz| < TT, as in the original model of Ref. [2]. For positive we do have | argz| < vr. 
Thus, we can find the OPE expression from the hadronic one, and vice versa, obtain the 
hadronic expression from the OPE. In particular, we can obtain its imaginary part, which 
corresponds to p(t). Part of the original structure of p(s) will be lost in the way, as the 
OPE misses non-analytic terms. The associated error is exponential with an exponent that 
vanishes when nj ^ (leaving aside the fact that one expects the series to be asymptotic). 
Therefore, the general structure of the solution would be 

p(s) = PpertXs) + (power like) + (exp. suppr.) , 

where only the last term does not follow from the OPE. Unfortunately this result has been 
obtained for an specific model. In general, one cannot get in a model independent way that 

P{s) = Ppert.is) H , 

where the dots mean a contribution that decays to zero faster than any power of 1/lns. 

5 Conclusions 

We have studied the constraints that the OPE imposes on models for current-current cor- 
relators inspired in the large limit and Regge trajectories. We have considered the case 
of scalar/pseudoscalar correlators. Assuming a model for the mass spectrum consisting of a 
linear Regge behaviour plus corrections in 1/n, we have obtained the logarithmic behaviour 
in n of the decay constants within a systematic expansion in 1/n. We have accomplished 
this by matching the hadronic and OPE expressions of the Adler function. The inclusion of 
1/n corrections to the decay constants is compulsory if one aims at going beyond the leading 
partonic result, as the 1/n terms are needed to produce logarithms of in the Euclidean. 
We have performed a numerical analysis of our results, in which we have seen the importance 
of the resummation of the In n terms of the decay constants, specially for low n. Our results 
show that it is possible to have a different slope, Bs/p, for the Regge behaviour of the scalar 
and pseudoscalar spectrum, and yet comply with all the constraints imposed by the OPE 
(including perturbative corrections). 

In order to check our setup in a controlled environment, where the results obtained with 
the OPE can be tested with the exact results, we have studied the preasymptotic effects of 
the scalar/pseudoscalar correlators in the 't Hooft model (QCD in 2 dimensions in the large 
Nc limit). This has allowed us to compute the 1/n corrections to the decay constants in 
the 't Hooft model for the first time. Actually, the connection between the OPE and the 
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decay constants provides us with a relatively easy way of finding 1/n corrections to hadronic 
matrix elements in the 't Hooft model. A direct analytic computation from the 't Hooft 
equation appears to be quite involved, however we have been able to confirm our results 
with the numerical evaluation from the 't Hooft equation. 

Finally, we have also considered a model at finite Nc t2j and modified it to allow for the 
inclusion of perturbative (logarithmic) corrections. We have then obtained the associated 
decay constants that are consistent with perturbation theory, though they are no longer 
trivially related with the perturbative expression of the imaginary part of the correlator. We 
have seen that consistency of this model with the OPE can be obtained with different slopes 
for the Regge trajectory associated to each channel. 
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A Appendix 

In this appendix we will present general, model- independent, matching formulas between 
the OPE and the hadronic calculation both for vector /axial- vector (for the notation in this 
case, see Ref. [TB]) and for scalar /pseudo-scalar correlators, within an expansion in powers 
of l/Ml 

The procedure is the same as in the main body of the text. Take the Adler function {Ax 
in the vector /axial- vector case, Bx in the scalar/pseudo-scalar case), transform it using the 
Euler-Maclaurin formula, and focus on the piece that can produce logarithms, namely, the 
integral. 



A.l Vector /Axial- vector correlators 

In this case we work with 

Flin) 



dn 



n* 



dM„ 



dM?, 



dn 



(79) 



Now, instead of assuming a model for we will just expand 



dM2 



dn 



Fl{n) 



Fx,o{Ml) + 



Fx,i{Ml) , Fx,2{M- 



M2 



+ 



Mi 



+ . . 



^0) 



Within this expansion, retracing the steps of section 3 in Ref. [12] the matching conditions 
between OPE and hadronic calculations are trivially obtained: 

. LO: FoAMl) = ^Imn^^-(M2) 
• NLO: Fi,x(M2) = 
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A. 2 Scalar/Pseudo-scalar correlators 

Here we work with 



dn 



^1) 



As the dimensions of these correlators are different from those of the vector/axial- vector 
correlators, we will use a different expansion, 

' FUn: ,) . Ml (f„(MJ; ,) + + +...). (82) 



(in 



?i n 



And the matching results are 

. LO: F,,x{Ml■^^) = -^hmW^iMl^i) 
• NLO:Fi,x(M2;;x) = 

. NNLO: FM,) ^ ^^,,,(^)^%,«.(M„^)M^^ 
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